Long-Time Effects in a Simulation Model of Sputter Erosion 
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A simple (2+1) dimensional discrete model is introduced to study the evolution of solid surface 
morphologies during ion beam sputtering. The model is based on the same assumptions about 
the erosion process as the existing analytic theories. Due to its simple structure, simulations of 
the model can be performed on time scales, where effects beyond the linearized theory become 
important. Whereas for short times we observe the formation of ripple structures in accordance 
with the linearized theory, we find a roughening surface for intermediate times. The long time 
behavior of the model strongly depends on the surface relaxation mechanism. 
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Introduction During the last years, two features 
of surface morphologies created by ion beam sputter- 
ing attracted particular attention: ripple structures 
on sub-micrometer length scales and self-afflne, rough 
surfaces Q. 

The formation of periodic ripple structures has been 
observed experimentally in amorphous materials [0, 
metallic crystals [pi 0] and semiconductors amorphized 
by the ion beam 1, |. Ripples are typically oriented 
perpendicular to the projection of the ion beam in the 
surface plane for small angles of incidence 6 (relative to 
the surface normal), whereas for larger angles 0, the ob- 
served ripple pattern is rotated by 90°. Surfaces eroded 
by ion bombardment may also exhibit self-afflne prop- 
erties J7J. With increasing ion fluence, a crossover from 
ripple structures to self-afflne, rough surfaces has recently 
been observed experimentally j|. 

Our present understanding of these features is based 
upon the work of Bradley and Harper (BH) H, who 
found that Sigmund's sputtering theory |J implies a cur- 
vature dependence of the sputtering yield. Based on BH a 
continuum theory of surface evolution by sputter erosion 
was formulated as an anisotropic Kuramoto-Sivashinsky 
(KS) equation @ with additive noise 0, H, JL3|. How- 
ever, there are strong indications from experiment [|[ [lj] 
that surface relaxation processes, are also important dur- 
ing pattern formation. Such processes have not yet been 
adequately included in existing theories. Our results will 
show that the long-time behavior of patterns depends 
crucially on details of surface relaxation. 

To investigate the analytic theory beyond the lin- 
earized regime, numerical integrations of the KS equa- 
tion have been performed |U| [|l6[ , which uncovered two 
markedly different long-time regimes, depending on the 
signs of the non-linear couplings. 

Computer simulations may be helpful in clarifying 
both the role of surface relaxation and of non-linear ef- 
fects. Two types of simulations have been performed 



up to now. Koponen et al jl7[ [l8| calculated collision 
cascades emerging from single ion impact within the bi- 
nary collision approximation. They find ripples in accor- 
dance with linearized BH theory JL8| , which appear both 
with and without additional surface relaxation processes. 
This indicates the presence of an ion-induced surface dif- 
fusion mechanism, which has also been predicted from 
BH fl3|| . The simulations did not yet reach time scales 
where the non-linear effects of the continuum theories 
could be analyzed. On the other hand, scaling proper- 
ties of the roughness of ion-irradiated surfaces have also 
been investigated within this approach (l7| . In a different 
approach, Cuerno et al [|l9| proposed a simple, discrete 
stochastic model with an update rule, which incorporates 
the ©-dependence of the sputtering yield and a simple 
curvature dependence of the erosion probability ad hoc. 
Within this model, it is possible to study the crossover 
from ripples to rough surfaces during the evolution of an 
irradiated 1-dimensional system. 

In the present letter, we introduce a simple, yet atom- 
istic Monte Carlo model on a lattice, which includes the 
same assumptions on energy deposition from ion impact 
as BH and we study the evolution of surface morpholo- 
gies beyond the linearized regime for two different types 
of surface relaxation mechanisms. 

Model We model energy deposition by directly im- 
plementing the result of Sigmund's sputtering theory ||, 
upon which BH and the existing analytical theories are 
based. The surface is described as a two-dimensional field 
of discrete time-dependent height variables h(x,y,t) on a 
square lattice of size Lx L. t denotes the time, measured 
in terms of ion fluence. In our results, we also indicate 
the corresponding amount of eroded material in terms of 
eroded monolayers (ML). Periodic boundary conditions 
are assumed, i.e. h(x + L, y, t) = h{x, y, t) = h(x, y + L, t). 
Initially the surface is flat, i.e. h(x,y,0) = ho. Two pro- 
cesses take place: erosion and surface relaxation. For 
each erosion step an ion is started at a random position 



in [0,L] x [0,L] above the surface and moved towards 
the surface with incidence angle 6 and an angle <\> with 
respect to the x-axis. After the ion has penetrated a 
distance a under the surface it stops and distributes its 
energy. A particle at the surface obtains the energy given 
by Sigmund's theory M : 



E(x',y',z') 



(2ir) 3 / 2 aiJ : 



■ exp 



0' + a) 2 x' 2 + y' : 



2a 2 
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(1) 

where (x',y',z') is the position of a particle in a local 
cartesian coordinate system of the ion, where the z' axis 
coincides with the ion trajectory. A particle at point r = 
(x, y, h) on the surface is removed, i.e. the height variable 
decreased by one, with a probability proportional to E(r) 
(see Fig. |). 
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FIG. 1: The model consists of a square field of discrete height 
variables h(x,y), corresponding to piles of h(x,y) particles at 
position (x,y). Left: Each ion impact is modeled by an dis- 
tribution describing the energy deposited by the ion. Atoms 
on the surface are removed with a probability proportional 
to the energy. Right: Surface diffusion, by decreasing height 
differences the energy is decreased. 

For surface relaxation, we use either activated surface 
diffusion (ASD) or Wolf- Villain (WV) type irreversible 
relaxation as simplified models of high-temperature and 
low-temperature relaxation mechanisms. ASD is imple- 
mented as a nearest neighbor MC hopping process with 
an energy E = (J/2) £<;, i? >(7*« - KJ)Y as in Ref. |j 
(see also right part of Fig. |l|) and WV relaxation is im- 
plemented according to Ref. |2j. For ASD a MC step is 
performed every 0.001L 2 erosion steps with temperature 
kfrT/J = 0.2. The irreversible WV process is completed 
for all surface atoms after every ion impact. System sizes 
up to L = 256 are considered. To avoid anisotropics, 
which may be induced by the square lattice, we have 
chosen the azimuthal angle <\> = 22°. The angle was 
varied between 0° and 80°. 

We use parameters a = 5.4, a = 3.3 and fi — 1.7 and 
e = (27r) 3 / 2 cr/i 2 (which makes the prefactor in Eq. (|l|) 
equal to one). One should keep in mind that our choice 
of e leads to rather high sputtering yields Y(Y « 7.0) 
compared to experiments like 0, where Y = 0.3 . . .0.5. 
According to BH, the ripple wavelength A scales like A ex 
Y~ x l 2 so that in our simulations we expect patterns with 
correspondingly smaller length scales. 



Results The model perfectly reproduces the 6 de- 
pendence of the sputtering yield of BH J22J and thus 
should not be used for too large 0. Wavelength ratios 
obtained from simulations ( A(45°)/A(30°) = 0.89(10), 
A(60°)/A(30°) = 0.84(10), A(70°)/A(30°) = 0.89(10) 
agree reasonably well with linear BH theory and with 
experiment [[|. 

Next, we present results on the evolution of the mor- 
phology of a surface relaxed by ASD to discuss the 
crossover between different time regimes. Afterwards we 
will show that the morphologies reached with WV relax- 
ation are completely different. 

For ASD, real space pictures of the height profiles of 
a 256 x 256 surface irradiated with ions under = 50° 
are shown in Fig. || for increasing ion fluences. Fig. || 
shows structure factors 5(k) = |/^(k)| 2 for wavevectors fey 
parallel and k± perpendicular to the projection of the di- 
rection of the ion beam onto the fiat surface. Here /i(k) is 
the Fourier transform of height fluctuations (h(x, y) —h) 
(h = ^2 X h(x,y)/L 2 being the average height). First 
signs of a structure appear after 1 ML of eroded mate- 
rial. For short times, after three MLs have been sputtered 
(t = 2.7 x 10 4 ions), a pattern with a typical length scale 
of 7 lattice spacings can be seen in the real space pictures. 
After 30 MLs have been eroded (t = 2.5 x 10 5 ) ripples ori- 
ented perpendicular to the direction of the ion beam are 
visible and maxima in the corresponding structure factor 
appear. For intermediate times (t = 2.6 x 10 6 correspond- 
ing to 300 MLs), the ripples start to disappear again. 
After long times (t = 9.3 x 10 7 or 10 4 MLs) the morphol- 
ogy approaches a rough surface with anisotropic statistics 
of height fluctuations. In the structure factor, the peak 
indicating the ripples has strongly decreased and S(k\\) 



k\ 



-2.83 



approaches a scaling behavior, S(k\\) 
lar crossover from ripples to a rough surface is also seen 
in experimentally obtained structure factors H. Please 
note that we also have tested that when switching of 
the sputtering mechanism, the structure factor obtained 
from pure surface diffusion follows at long times the usual 
S(k)~ k~ 2 behavior §. 

The different time regimes are seen most clearly 
in the time evolution of the roughness W 2 (t) = 
l/L 2 ^ xy (h(x,y,t)-h(t)) 2 (see Fig. |). For short times 
the data exhibit a positive curvature. During that period 
the ripple structure is formed. The crossover to interme- 
diate times appears rather abruptly in accordance with 
numerical solutions of the noisy KS equation ]16[ . A fit 
with an algebraic growth law W ~ t@ for intermediate 
times reveals /3 = 0.168(2). For longest times, a strong 
increase in W is observed. We assume that the end of 
the algebraic regime is marked by finite size effects, which 
remain to be studied in detail. 

For WV relaxation and for short times, we again find 
ripples which behave in accordance with linear BH the- 
ory. The long-time morphologies, however, are com- 
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FIG. 2: Surface of a 256 x 256 system after t ions have hit 
the system. G = 50°. Upper row left: t — 2.7 x 10 4 (corre- 
sponding to 3 MLs of particles been eroded), upper row right: 
t = 2.5 x 10 5 (30 MLs), lower row left: t = 2.6 X 10 6 (300 
MLs), lower row right: t = 9.3 x 10 7 (10 4 MLs). The bars 
indicate the direction of the ion beam. 



pletely different from ASD relaxed surfaces as can be 
seen in Fig. |5|. Generally, WV produces stable patterns 
of pronounced ripples up to the longest times simulated. 
For © = 50° the ripples are tilted by « 17° against the di- 
rection perpendicular to incidence, whereas for 6 = 70° 
striplike domains oriented approximately parallel to in- 
cidence appear, which contain ripples with the same tilt 
angle of ~ 17°. 

In conclusion, we have presented a simple discrete 
model for (2+l)-dimensional surfaces subjected to ion- 
beam sputtering. The model is based on assumptions, 
which also underly BH theory and continuum theories. 
The simplicity of the model allows us to simulate rather 
large systems over time scales, where non-linear effects 
become apparent. We have used the model to study dis- 
tinct time regimes of surface evolution under the influ- 
ence of two different surface relaxation mechanisms (ASD 
and WV). For short times, we always find ripple forma- 
tion in accordance with linear BH theory. 

For longer times, the surface morphologies emerging 
from ASD and from WV become qualitatively different. 
For ASD, the ripple wavelength increases with increasing 
ion fiuence, and scaling behavior in the structure factor 
appears. Ripples get blurred and the surface morphology 
becomes rough, but stays anisotropic. Replacing ASD by 
WV surface relaxation, we find pronounced ripples in the 
long-time regime. Please note that in the earlier detailed 
simulations of the sputtering process [|1|, basically the 
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FIG. 3: Structure factor S(k\\) (upper panel) and S(k±) 
(lower panel) for different ion fluences and for G = 50°. The 
solid line in the upper panel represents a k~ 2 ' 83 behavior. The 
data for the largest ion fiuence is an average over 10 indepen- 
dent runs, while the other data points are from 300 runs. 



same result was obtained even when no surface diffu- 
sion was present. It was assumed that this is due to the 
smoothing mechanism introduced by the detailed mod- 
eling of the sputtering itself. In our simple model, we 
need an additional surface relaxation process to observe 
ripples, but we have found for short times no significant 
differences between the WV and ASD cases. Hence, the 
non-linear long time regime allows a much better com- 
parison of different surface diffusion mechanisms. 

WV relaxation and ASD may be considered as the 
simplest type of relaxation mechanisms for low temper- 
atures and for high temperatures, respectively. Please 
note that the aim of this work is not the evaluate 
the relaxation mechanisms itself, but to introduce a 
simple model which allows to study long-time effects 



10° 



i<r 



.no 1 



10 u 











_.....---— 


v^" v 


r 










- 




,/io 4 












1 




1 








10 3 

10 2 


. 0=80° 

©=30° 

©=0° 


// 


•'TV*-* 








10 1 
10° 

1( 


- ^\^' 










f 10 5 




10 7 







10" 



10° 



10° 
t/lons 



10' 



10° 



FIG. 4: Surface roughness W 2 as a function of number t of 
ions for G = 50°. 




FIG. 5: Surfaces after long time sputtering with Wolf- Villain 
type surface diffusion, left: G = 50° right: G = 70° 



and their dependence of the surface diffusion mecha- 
nism. Other, more specific surface relaxation mecha- 
nisms ||, H H ||, |7|, H, |9[ |0| involving diffusion 
barriers may be studied within our model using appro- 
priate kinetic MC steps for surface relaxation, especially 
the model can be easily extended to describe different 
crystal structures and include surface anisotropics. 
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